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Wu  (1985)  and  Zeng  et  al.  (1991,  1993)  have  shown  that  the  energy  transfer  theory 
can  successfully  describe  the  seismic  body  waves  scattering  processes  in  a  random 
scattering  and  absorption  full-space  medium.  In  this  study,  we  shall  first  use  the 
energy  transfer  equations  to  study  the  seismic  surface  waves  scattering  processes 
and  the  surface  wave  with  body  waves  scattering  processes  in  order  to  Interpret  the 
Su  and  Akl's  observed  results.  We  shall  then  extend  this  method  to  study  the 
explosions  and  earthquakes  at  regional  distance. 
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Study  on  the  Discrimination  of  Explosions  and  Earthquakes 

AT  REGIONAL  DISTANCE  BY  USING  CODAQc*  METHOD 

Xiaofei  Chen  and  Keiiti  Aki 

Department  of  Geological  Sciences,  University  of 
Southern  California,  Los  Angeles,  CA  90089-740 

Contract  Number.  F49620-93-1-0016 


Objective 

The  objective  of  this  study  is  to  discriminate  the  explosions  and  earthquakes  at  regional 
distance  by  using  the  coda  wave  method.  In  an  earlier  study,  Su  and  Aki  (1991)  found  a 
significant  difference  in  coda  attenuation,  Qc'^  for  quarry  blasts  and  earthquakes  at  frequencies 
of  1.5  and  3  Hz  for  lapse  time  less  than  30  seconds  (see  Fig.  1),  and  they  suggested  that  such  a 
significant  difference  in  Qc*'  may  be  attributed  to  the  seismic  surface  wave's  contributions.  To 
interpret  the  observed  seismic  coda  waves  we  need  to  consider  the  seismic  surface  wave 
scattering  processes.  Wu  (1985)  and  Zeng  et  al.  (1991,  1993)  have  shown  that  the  energy 
transfer  theory  can  successfully  describe  the  seismic  body  waves  scattering  processes  in  a 
random  scattering  and  absorption  full-space  medium.  In  this  study,  we  shall  first  use  the  energy 
transfer  equations  to  study  the  seismic  surface  waves  scattering  processes  and  the  surface  wave 
with  body  waves  scattering  processes  in  order  to  interpret  the  Su  and  Aid's  observed  results.  We 
shall  then  extend  this  method  to  study  the  explosions  and  earthquakes  at  regional  distance. 


Research  accomplished 

1.  Energy  transfer  theory  for  Rayleigh  wave 

In  this  section,  we  shall  consider  a  simple  case  in  which  the  background  medium  is  a 
homogeneous  half-space.  The  only  surface  wave  in  this  case  is  the  Rayleigh  wave.  In  this 
section,  we  neglect  the  conversions  between  the  body  waves  and  the  Rayleigh  wave,  and 
consider  only  energy  distribution  of  Rayleigh  wave  in  an  absorptive  and  scattering  medium. 

The  energy  density  at  x,  and  co,  due  to  a  point  steady  source  at  x,  can  be  written  as 

E'^(x,to)  =  G(x,xs)expr-ii*  lr-rjl]£g(Xs,(i))  +  ^  a(Xj)G(x,Xs)exp[-'n'^  lr-rjl]E^(Xj.ci)) ,  (1) 

j 

where,r=  x-ez(x*ez)  is  a  horizontal  vector,  and  tj"  is  the  anenuation  coefficient  that  includes  the 

intrinsic  (Tif )  and  scattering  attenuation  (1)5 ),  i.e.,  "n"  =  Tif  +  ti5.  In  eq.(l),  E^x,Q))  is  the  seismic 

energy  per  unit  volume  carried  by  Rayleigh  wave  at  point  x  and  frequency  co.  The  first  term  in 
the  right-hand-side  of  (1)  represents  the  direct  Rayleigh  wave  energy  radiated  fiom  source,  and 
the  second  term  in  the  right-hand-side  of  (1)  describes  the  total  scattering  energy  from  all 

scatterers  (xj).  Where  o  is  the  scattering  cross-section,  G(x,Xj)  and  exp[-T]"lr-rjl]  describes  the 
geometrical  spreading  and  attenuation,  respectively.  It  should  be  not^  that  the  attenuation  is 
caused  only  by  the  horizontal  propagation,  since  we  only  consider  surface  wave  in  this  section. 


For  weak  scattering  cases,  the  scattering  cross-section  a  can  be  analytically  derived  by  using 
Bom  approximation  (Aki  and  Richards,  1980).  Here,  for  simplicity,  we  assume  that  o  is  a  con¬ 
stant.  The  geometrical  spreading  funcdon  G(x,x'),  however,  can  be  determined  by  the  consid¬ 
eration  of  energy  conservation. 

Geometrical  Spreading  Function  of  Rayleigh  Wave: 

The  direct  Rayleigh  wave  (Rayleigh  wave  in  the  background  medium)  has  the  following 
form  (Aki  and  Richards,  1980), 

U'*(x,co)  =  exp{i[kRlr-rsl+Jt/4] ).  (2) 

8(cr)^Ii  V  (jt/2)kRlr-rsl 

Where, 

Pr(x)  =  ri(2,©)  Cr  +  ir2(z,co)ez  ,  (2a) 


with  Cr  =  (r-rs)/lr-rsl.  The  corresponding  energy  density  is 

ES(x,0))  =  i^0)2|U“(x,(0)|2  = - ^ IPr(Xs).F(Xs)|2  .  (3) 

2  [8(cR)2li]27ckR  lr-r,l 

For  steady  state,  the  energy  rate  should  be  a  constant  for  any  closed-surface  that  encloses  the 
source,  namely. 


{E2(x,0))vr}  *8  dS  =  constant. 


for  Xs  inside  S. 


(4) 


Where,  s  is  the  normal  vector  of  the  surface  element  dS,  and  vr  is  the  velocity  vector  of  the 
Rayleigh  wave.  Since  the  surface  wave  propagates  along  the  horizontal  direction,  we  consider  a 
cylindrical  surface  that  leads  to  a  simple  surface  integration. 


{ES(x,©)vR)«sdS 


RdO  e5(x,(i))vr  . 


Where,  R  =  Ir-fjl.  Substituting  (3)  in  (4),  we  find 

e5(x,)=^| 

=  — - IF(Xs)|2{iri(Zs)nr|2-»-(r2(zs)n/)  .  (5) 

8(cR)4likR 


{E?(x,(d)vr)  •sdS 


2(o^p  f 


IPR(z)|2dz  IPr(Xs)*F(Xs)I^ 


Using  the  above  result,  eq.  (3)  can  be  rewritten  as 


(6) 


Where. 


gR(Z,(0)  =  — 


IPr(z)|2 


I  IPR(z-)|2dz*  I 


[ri(z.co)^  +  r2(z.Q))^] 

» 

[ri(z'.(o)^  +  r2(z‘,(i))^]dz 


(6a) 


Thus,  we  obtain  the  geometrical  spreading  function  of  Rayleigh  wave  energy  propagation  as 

G(x.x)  (7) 

27clr-rl 

From  equation(6a),  we  can  verify  the  following  identity, 

gR(z,0))dz  =  1  .  (7a) 


I 


Integral  equation  of  scattering  Rayleigh  wave  in  random  medium: 

We  are  now  ready  to  set  up  the  basic  integral  equation  for  Rayleigh  wave  scattering  in 
random  medium.  Substituting  eq.(7)  in  (1),  and  assuming  that  the  random  scatterers  are 
uniformly  distributed  and  can  be  described  by  a  continuous  distribution  no(density  of  scatterers), 
we  obtain  the  following  integral  equation. 


E^x.O))  =  e?(x.a))  +  (  Titi5(^:^e-^“'r-f''E*(x,©)dV(x’)  .  (8) 

2nlr-rsl  J  27tlr-r’l 

Where,  =  noO.  Eq.(8)  is  similar  to  the  energy  transfer  equation  for  body  wave  (Wu,  1985; 
Zeng  et  al.,  1991),  but  with  different  geometrical  spreading  functions  and  the  attenuation  factors 
that  are  due  to  the  characterizations  of  surface  wave  propagation.  The  dependence  on  Ir-r'l 
indicates  the  horizontal  propagation  property  of  the  surface  wave,  whereas,  the  gR(z,co)  expresses 
the  depth  distribution  of  surface  wave  energy.  Integral  equation  (8)  can  be  solved  by  using 
spatial  domain  Fourier  transform  method.  First,  we  note  that  equation  (8)  can  be  rewritten  as 

E*(x,©)  =  gR(z,(o)E  j*(r,0)) .  (9) 

Inserting  equation  (9)  into  (8)  and  using  the  identity  (7a),  we  obtain 

E5(r,co)  =  e5(Xs,a))  +  (  r]f  Ej(r',o))  dZ(r’) . 

2rar-rsl  2n\r-r{ 


(10) 


Where,  2(r)  =  {(x,y)l-«xx<+oo,  -oo<y<+oo)  .  The  corresponding  Fourier  transform  is 


E2(k.(D)  =  G2(k)eS(x,.(0)  +  Tis'*G5(k)E5(k.o}) . 


(11) 


Where,  we  assumed  r*  =  0,  and  G2(k,0))  is  given  by 


G2(k,co)  =11  dxdy  exp(ik»r)  =  I  Jo(kr)e-’i"^dr  =  1 

JJ  27tr  I  Vk2  +  (n*)2 


(11a) 


Substituting  this  result  in  eq.(ll),  then  taking  the  inverse  Fourier  transform  over  k,  we  finally 
obtain 


E'‘(x,a))  =  gR(z,co)PR(r,Ti*,Ti*)  eS(Xs,o>) . 

Where  PR(r,Ti*^,T\s*)  is  the  inverse  Fourier  transform  of  G2(k,a)) , 


(12) 


PR(r,'n’‘,Tl?)  =  ^ 


Jn(kr)— —  kdk  - 

Vk2  +  -  Ilf 


(12a) 


Solution  (12)  has  a  clear  physical  meaning.  The  term  e5(Xs,(o),  as  defined  earlier,  contains  the 
seismic  source  information.  Function  gR(z,0))  represents  the  depth  dependence  of  Rayleigh  wave 
energy,  whereas  the  function  PR(r,Ti‘‘,TiJ)  describes  the  propagation  and  attenuation  processes  of 
the  Rayleigh  wave  energy.  PR(r,il‘‘,Tis)  can  be  obtained  by  evaluating  integral  (12a).  For  a  pure 

absorption  medium  (Tls=0),  we  find  PR(r,Ti‘*,0)  =  e*hV27cr .  This  is  consistent  with  our  direct 
Rayleigh  wave  (solution  in  background  me^um).  For  general  absorption  and  scattering  media, 
we  can  numerically  evaluate  the  propagation  and  attenuation  function  PR(r,Ti",T|,"). 


It  is  noted  that  equation  (8)  describes  a  stationary  energy  transfer  process  ,  i.e.,  E**(x,©) 
represents  the  amplitude  of  the  spectrum  of  scattered  Rayleigh  wave.  To  obtain  the  time-history 
of  the  scattered  Rayleigh  energy  for  a  given  frequency  w  and  the  observational  point  x,  we 

introduce  a  time-delay  phase  factor  of  “exp(-iQlr-r’l/cR)”,  then  take  inverse  Fourier  transform 

over  frequency  Q.  During  the  inverse  Fourier  transformation,  the  frequency  to  is  kept  as  a 
constant  which  is  defined  as  center  frequency.  Thus  we  can  obtain  the  time-history  of  scattered 

energy  for  a  given  center  frequency  to.  Fig.  2  shows  the  energy  distributions  for  the  case  of 

Tif=0.01,  q^=0.02  and  a)=lHz.  As  we  expected  the  scattered  energy  arrived  at  time  of  t=lr-rsl/cR. 


2.  Energy  transfer  theory  for  S  and  Rayleigh  waves’  scattering  processes 

In  the  preceding  section,  we  have  considered  only  the  scattering  processes  of  Rayleigh 
wave,  and  neglected  the  body  waves’  conversions.  In  this  section,  we  shil  consider  the  coupling 


effect  of  body  waves  with  surface  wave.  For  simplicity,  we  consider  only  the  coupling  process 
between  the  S  wave  and  Rayleigh  wave.  The  energy  transfer  equation  described  such  an  coupled 
scattering  process  can  be  written  as  follows, 

E"(x,(d)  =  G*(x,Xs)e-^"**’-'‘Jeo  '‘(x,to)  ( 1 3a) 

+  j  G*(x,x‘)e-’i"''--fl {tIs^EV,©)  +  Ti,®E*(x',©)}dV(x) 

E^x,©)  =  G®  (x,Xs)e-^*'''-f*'eo®  (x,©)  ( 1 3b) 

+  1  G*(x.x)e-^*'*-*4ils“EV.co)+q,“EV.£o)}dV(x')  . 

Where, 

ES(x,©):  seismic  energy  carried  by  S  wave  per  unit  volume  at  x; 

;  total  attenuation  coefficient  for  Rayleigh  wave,  =  qf  +  qP  +qp; 
qP:  Rayleigh  to  Rayleigh  waves  scattering  coefficient; 
qf :  S  to  Rayleigh  waves  scattering  coefficient; 
q^  :  total  attenuation  coefficient  for  S  wave,  q^  =  qf  +  qP  +qf®; 
qf  :  absorption  coefficient  for  S  wave; 
qP  :  Rayleigh  to  S  waves  scattering  coefficient; 
qf  ^  ;  S  to  S  waves  scattering  coefficient; 

eoCxs.co) :  total  S  wave  energy  rate  radiated  from  source  divided  by  the  velocity  of  S  wave; 
G*^(x,x') ;  the  geometrical  spreading  function  for  Rayleigh  wave  given  by  equation  (7); 

G^(x,x'):  the  geometrical  spreading  function  for  S  wave  in  a  half-space  medium,  and  it  can 
be  approximated  by 


G*(x,x’ )  =  — 1 —  + - 1 - . 

47tlx-xl  4jtlx-(x’)*l 

Where,  (x’)*  is  the  image  point  of  x’  with  respect  to  the  free  surface  z=0,  i.e.,  (x’)*=(x’,  y’,-z’). 

Equations  (13a)  and  (13b)  can  be  further  simplified.  Making  a  two-dimensional  Fourier 
transform  over  the  horizontal  variable  (x,y),  equation  can  be  reduced  to. 


E*(K,z,©)=  G'*(K,z;z,;©)eg(x„©) 

+ 1  g'‘(K,z;z’;©)  {qPf‘(K,z’,©)  +  qs='*E®(K,z’,©))dz’  , 


(14a) 


E®(K,z,©)  =  G*(K,z;Zs;©)4(Xs,©) 


(14b) 


Where,  K  is  the  horizontal  wave-numbcr  and 


and 


*(K,z;z';Ci)  G(r,z;z';i 
j^(K,z;(0)=g(z,(0)j 


ci))exp(-Ti* V r2+(z-z’)^  )Jo(Kr)rdr, 


exp(-‘n’^r)Jo(Kr)dr . 


Solving  the  coupled  equations  (14a)  and  (14b),  we  can  obtain  the  solution  of  energy  transfer 
processes  of  S  and  Rayleigh  waves.  The  time-history  of  total  scattering  energy  can  bt 
determined  by  introducing  proper  time-delay  phase  factors  to  these  coupled  equations.  Figure. 
(3)  show  the  time  history  of  total  scattering  energy  of  S  and  Rayleigh  waves,  for  the  case  of 

Tif=0.01(km)’l,  TiP=0.01(km)‘l,  T)f^=0.02(km)‘l,  T)f=0.03(km)‘l,  ■nP=0.02(km)'l  and  Tif*^= 
0.01  (km)' 1  for  various  source  depths  and  center  frequencies.  Our  results  indicate  that  the 
shallower  the  source  depth  and  the  lower  the  frequency  is,  the  contribution  of  scattered  Rayleigh 
waves  is  larger.  This  is  consistent  with  the  observed  results  of  Su  and  Aid  (1991)  as  shown  in 
Fig.  I,  where  earthquakes  have  an  average  focal  depth  of  8  km,  whereas  the  average  focal  depth 
of  quarry  blasts  is  about  10m.  Therefore,  the  latter  contains  much  Rayleigh  waves’  contributions 
while  the  former  dominated  by  S  waves  only.  We  also  found  that  the  contribution  of  Rayleigh 
waves  decreases  as  the  center  frequency  increasing.  The  suppressing  of  surface  wave 
contribution,  however,  is  not  strong  enough  to  directly  fit  the  observed  results  shown  in  Fig.  1 . 
This  indicates  strong  attenuation  of  surface  waves  at  higher  frequencies  due  to  absorption. 


Conclusions  and  future  studies 

To  date,  we  have  developed  the  energy  transfer  theory  for  Rayleigh  and  S  wave 
propagation  in  a  random  scattering  and  attenuation  half-space  medium.  The  solutions  of  our 
energy  transfer  equations  can  qualitatively  explain  the  Su  and  Aki’s  observed  results  about  the 
local  earthquakes  and  quar^  blasts  (1991).  Our  study  indicates  that  the  surface  wave  scattering 
become  important  for  the  events  with  shallower  focal  depth,  for  instance,  the  quarry  blasts  and 
nuclear  explosions.  For  the  explosions  at  a  regional  distance,  we  expect  similar  conclusions.  To 
solve  such  scattering  problem  at  a  regional  distance  we  shall  extend  our  energy  transfer  equations 

to  more  general  cases.  At  the  same  time  we  shall  calculate  the  coda  of  the  Chinese  nuclear 
explosion  and  the  nearby  earthquakes  by  using  the  CDSN  data.  We  expect  that  a  similar 
discrimination  between  explosions  and  earthquakes  as  observed  by  Su  and  Mi  for  local  events 
(1991)  can  be  found  in  the  events  at  regional  distance,  and  our  energy  transfer  theory  can  offer  a 
physical  basis. 


References 


Aki,  K.  and  P.  G.  Richards  (1980).  Quantitative  Seismology:  Theory  and  Methods,  W.  H. 
Freeman,  San  Francisco. 

Su,  F. ,  K.  Aki,  and  N.  N.  Biswas  (1991).  Discriminating  quarry  blasts  from  earthquakes  using 
coda  waves.  Bull.  Seism.  Soc.  Am.,  81, 161-178. 

Wu,  R.  S.  (1985).  Multiple  scattering  and  energy  transfer  of  seismic  waves—separation  of 

scattering  effect  from  intrinsic  attenuation,  I.  Theoretical  modeling,  Geophys.  J.  R.  Astr. 
Soc.,  82,  57-80. 

Zeng,  Y.  H.,  F.  Su,  and  K.  Aki  (1991).  Scattering  wave  energy  propagation  in  a  random 
isotropic  scattering  medium,  1.  Theory,  J.  Geophys.  Res.,  96, 601-620. 

Zeng,  Y.  H.  (1993).  Theory  of  scattered  P  and  S  waves  energy  in  a  random  isotropic  scattering 
medium,  BSSA,  83, 1264-1276. 


(1/Q)x1000 


Figure  1.  Coda  Q*^  vs.  lapse  time  obtained  using  the  local  earthquakes  and 
quarry  blasts  in  the  south  central  Mojave  Desert  area,  each  open  circle  on  the  plot 
represents  tme  measurement  fen:  a  particular  seismogram  on  a  time  window  of  34  sec 
for  fimquency  1.5  Hz,  25Hz  for  3  and  20sec  fenr  6Hz  and  12Hz.  The  solid  line 
connects  the  mean  points  (solid  circles)  calculated  by  averaging  the  individual 
measurements  in  each  8  sec  time  interval  with  4  sec  overiapping  at  the  adjacent  mean 
points.  The  standard  error  of  the  mean  is  also  shown  for  each  mean  point  (from  Su 
and  Aid,  1991). 
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Figure  2.  Scattered  Rayleigh  energy  for  the  case  of  Ttf  ^.03(km)‘* , 
Tif =0.02(kin)*^  and  <o=lHz.  Where  r=30km  and  CKs2.8kiii/sce. 


Total  Scattered  EfWBy(R>S)  Total  Scattered  Enersy  (R^) 


zs=10ktn  (dashed  line),  (b)  For  the  cases  of  a)=I.5Hz  (solid  line) .  O)=3.0Hz  (dashed  line)  and  o)=6Hz 
(dotted  line). 


